Cosmological limit on the neutrino mass 
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We have performed a careful analysis of constraints on the neutrino mass from current cosmological 
data. Combining data from the cosmic microwave background and the 2dF galaxy survey yields an 
upper limit on the sum of the three neutrino mass eigenstates of ^ m v < 3 eV (95% conf.), without 
including additional priors. Including data from SNIa observations, Big Bang nucleosynthesis, and 
HST Hubble key project data on Ho tightens the limit to < 2.5 eV (95% conf.). We also 

perform a Fisher matrix analysis which illustrates the cosmological parameter degeneracies affecting 
the determination of ^ m„. 

PACS numbers: 14.60.Pq,95.35.+d,98.80.-k 



I. INTRODUCTION 

The absolute value of neutrino masses are very dif- 
ficult to measure experimentally. On the other hand, 
mass differences between neutrino mass eigenstates, 
(mi, 77i2, ^73), can be measured in neutrino oscillation ex- 
periments. Observations of atmospheric neutrinos sug- 
gest a squared mass difference of 5m? ~ 3 x 10 -3 eV 2 
jl], D . While there are still several viable solutions to the 
solar neutrino problem the so-called large mixing angle 
solution gives by far the best fit with 8m 2 ~ 5 x 1CP 5 
eV 2 |,@. 

In the simplest case where neutrino masses are hi- 
erarchical these results suggest that mi ~ 0, 777,2 ~ 
5m so iar, and 777,3 ~ ^atmospheric- If the hierarchy is 
inverted || ^, [?J |8[ M, [l(]] one instead finds 7773 ~ 0, 

777 2 ~ ^atmospheric, and 777i ~ 5m atro o S pheric- HoW- 

ever, it is als o p ossible that neutrino masses are de- 
generate @ || g 0, || £| || |§ || |lj, 

777i ~ w>2 ~ "73 3> ^^atmospheric, in which case oscillation 
experiments are not useful for determining the absolute 
mass scale. 

Experiments which rely on kinematical effects of the 
neutrino mass offer the strongest probe of this overall 
mass scale. Tritium decay measurements have been able 
to put an upper limit on the electron neutrino mass of 
2.2 eV (95% conf.) [p2[ . However, cosmology at present 
yields an even stronger limit which is also based on the 
kinematics of neutrino mass. Neutrinos decouple at a 
temperature of 1-2 MeV in the early universe, shortly be- 
fore electron-positron annihilation. Therefore their tem- 
perature is lower than the photon temperature by a fac- 
tor (4/11) 1 / 3 . This again means that the total neutrino 
number density is related to the photon number density 
by 
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(1) 



Massive neutrinos with masses 777 ;§> To ~ 2.4 x 10 



eV are non-relativistic at present and therefore con- 
tribute to the cosmological matter density p3|, |24|, |25| 



92.5 eV 



(2) 
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calculated for a present day photon temperature Tq = 
2.728K. Here, ^ m v = m\ +7772 + 7773. However, because 
they are so light these neutrinos free stream on a scale of 
roughly k ~ 0.03777 o v^m /2 h Mpc -1 g § §. Below 
this scale neutrino perturbations are completely erased 
and therefore the matter power spectrum is suppressed, 
roughly by AP/P ~ -80 !/ /n ro f§. 

This power spectrum suppression allows for a deter- 
mination of the neutrino mass from measurements of the 
matter power spectrum on large scales. This matter spec- 
trum is related to the galaxy correlation spectrum mea- 
sured in large scale structure (LSS) surveys via the bias 
parameter, b 2 = P g (k)/P m (k). Such analyses have been 
performed several times before J29|, ||(|, most recently 
using data from the 2dF galaxy survey J3l[ . This inves- 
tigation finds an upper limit of 1.8-2.2 eV for the sum 
of neutrino masses. However, this result is based on a 
relatively limited cosmological parameter space. 

It should also be noted that, although massive neu- 
trinos have little impact on the cosmic microwave back- 
ground (CMB) power spectrum, it is still necessary to 
include CMB data in any analysis in order to determine 
other cosmological parameters. 

In the present paper we perform an extensive analysis, 
carefully discussing the issue of parameter degeneracies. 
The next section is devoted to a Fisher matrix analy- 
sis of the problem which establishes possible parameter 
degeneracies and yields a general idea of the precision 
with which the neutrino mass can be measured. Section 
III describes a full numerical likelihood analysis of data 
from CMB and LSS which yields a robust limit on the 
neutrino mass. Finally, section IV contains a discussion 
and conclusion. 
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II. FISHER MATRIX ANALYSIS 

Measuring neutrino masses from cosmological data is 
quite involved since for both CMB and LSS the power 
spectra depend on a plethora of different parameters in 
addition to the neutrino mass. Furthermore, since the 
CMB and matter power spectra depend on many differ- 
ent parameters one might worry that an analysis which 
is too restricted in parameter space could give spuriously 
strong limits on a given parameter. Therefore, it is desir- 
able to study possible parameter degeneracies in a simple 
way before embarking on a full numerical likelihood anal- 
ysis. 

It is possible to estimate the precision with which the 
cosmological model parameters can be extracted from a 
given hypothetical data set. The starting point for any 
parameter extraction is the vector of data points, x. This 
can be in the form of the raw data, or in compressed form, 
typically the power spectrum (Ci for CMB and P(k) for 
LSS). 

Each data point has contributions from both signal and 
noise, x = ^signal + a; no i sc . If both signal and noise are 
Gaussian distributed it is possible to build a likelihood 
function from the measured data which has the following 
form |3l 



£(9) oc exp 



\x\C{Q\ 



(3) 



where = (CI, ftj,, H , n s , r, . . .) is a vector describing 
the given point in model parameter space and C(9) = 
(xx T ) is the data covariance matrix. In the following we 
shall always work with data in the form of a set of power 
spectrum coefficients, Xi, which can be either C\ or P(k). 

If the data points are uncorrelated so that the data 
covariance matrix is diagonal, the likelihood function can 
be reduced to C oc e~ x > , where 



N n 



E 



(*^i,obs theory ) 

a(xi) 2 



(4) 



is a x 2 -statistics and N max is the number of power spec- 
trum data points [32] ] . 

The maximum likelihood is an unbiased estimator, 
which means that 



<e> = e„. 



(5) 



Here 0o indicates the true parameter vector of the under- 
lying cosmological model and (9) is the average estimate 
of parameters from maximizing the likelihood function. 

The likelihood function should thus peak at 9 ~ 9o, 
and we can expand it to second order around this value. 
The first order derivatives are zero, and the expression is 
thus 



+ >>-0) 



v.min ■ / j 

i,3 



dx k dx k 



dQi 89, 



(^/c,obs 



d 2 x k 



09,89, 



(6) 



r 



where i,j indicate elements in the parameter vector 9. 
The second term in the second derivative can be expected 
to be very small because (xk,obs — x k) is in essence just 
a random measurement error which should average out. 
The remaining term is usually referred to as the Fisher 
information matrix 



d 2 x 2 

89,89, 



Fi, 



1 8x k 8x k 



E a (x k )2~89~d9~ 



(7) 



The Fisher matrix is closely related to the precision with 
which the parameters, 9i, can be determined. If all free 
parameters are to be determined from the data alone 
without any priors then it follows from the Cramer- Rao 
inequality |33| that 



(8) 



for an optimal unbiased estimator, such as the maximum 
likelihood @. 

In order to estimate how degenerate parameter i is 
with another parameter, j, one can calculate how a(9i) 



changes if parameter j is kept fixed instead of free in the 
analysis. Starting from the 2x2 sub- matrix 



one then finds 



&j fixed (9i 



V (s- 1 ) 

We therefore define the quantity 

_ CTj fixed (^i) 



< 1 



(9) 



(10) 



(11) 



as a measure of the degeneracy between parameters i and 
j- 

In order to perform an actual calculation we use the 
most present data from CMB and LSS. 

CMB data set — Several data sets of high precision 
are now publicly available. In addition to the CGBE J3f| 
data for small I there are data from BOGMERANG [p6||, 
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MAXIMA (3J, DASI |3§ an d several other experiments 
p9| , |o). Wang, Tegmark and Zaldarriaga [|9) have com- 
piled a combined data set from all these available data, 
including calibration errors. In the present work we use 
this compiled data set, which is both easy to use and in- 
cludes all relevant present information. Altogether there 
are 24 CMB data points in this compilation. 

LSS data set — At present, by far the largest survey 
available is the 2dF of which about 147,000 galaxies 
have so far been analysed. Tegmark, Hamilton and Xu 
p2| have calculated a power spectrum, P(fc), from this 
data, which we use in the present work. The 2dF data 
extends to very small scales where there are large effects 
of non-linearity. Since we only calculate linear power 
spectra, we use (in accordance with standard procedure) 
only data on scales larger than k — 0.2ft, Mpc" 1 , where 
effects of non-linearity should be minimal. Making this 
cut reduces the number of power spectrum data points 
to 18. 

For calculating the theoretical CMB and matter power 
spectra we use the publicly available CMBFAST package 
P3fl . As the set of cosmological parameters we choose 
f2 m , the matter density, f2fc = 1 — Q m — fl\ — f2„, the 
curvature parameter, f^, the baryon density, Ho, the 
Hubble parameter, n s , the scalar spectral index of the 
primordial fluctuation spectrum, r, the optical depth to 
reionization, Q, the normalization of the CMB power 
spectrum, &, the bias parameter, and f2„, the neutrino 
density. In all cases we take the number of massive neu- 
trinos to be 3. The reason is that if neutrinos are to have 
an impact on CMB and matter spectra their masses must 
be much larger than the mass splitting inferred from at- 
mospheric neutrino observations (dm ~ 0.05 — 0.1 eV), 
and therefore neutrino masses will be degenerate with 

mi - m 2 ~ m 3 > (^atmospheric- 

In principle one might include even more parameters 
in the analysis, such as r, the tensor to scalar ratio 
of primordial fluctuations. However, r is most likely 
so close to zero that only future high precision exper- 
iments may be able to measure it. The same is true 
for other additional parameters. Deviations from the 
slow-roll prediction of a simple power-law initial spec- 
trum p4 145[ ftc} fl7| or additional relativistic energy den- 
sity g]g7g|jg ||, §4], |H Hi could be present. 
However, such effects only appear in cosmological models 
which are more complicated than the "standard" ACDM 
model. The parameters we use fully describe the features 
of the simplest working model. 

In the end one can check the consistency of the nu- 
merical parameter extraction by calculating the \ 2 P er 
degree of freedom. In Section III we find that the best fit 
is consistent with expectations and therefore it is unlikely 
that there are other parameters significantly affecting the 
power spectra. 

Fig. 1 shows the matrix r,j, calculated for the WTZ 
+ 2dF data set errors, around a reference cosmological 
model with parameters fl m = 0.3, fl\ = 0.7, fl},h 2 = 
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FIG. 1: Values of the parameter ry, defined in Eq. ( |ll| ) 



the ACDM concordance model. Note that r u is an ill- 
defined quantity. For plotting purposes we simply put 
rn = 1, but this has no physical significance. 

From this general matrix one can study the m„-related 
degeneracies more closely. Clearly the two parameters 
most degenerate with m v are f2 m , the matter density, 
and b, the bias. This is not surprising because a non- 
zero neutrino mass has little effect on the CMB acoustic 
peaks and a large effect on the matter power spectrum 
on scales below the free streaming scale. 

Therefore, any parameter which behaves in a similar 
way will cause degeneracy. This is indeed the case for Q m , 
if this parameter is changed while keeping flk fixed there 
is little effect on CMB. On the other hand, changing f2 m 
changes the normalization of the matter power spectrum 
at small scales relative to large scales. Changing b also 
mimics a non-zero neutrino mass. The reason is that 
present day LSS data have large error bars around the 
free-streaming scale for light neutrinos in the eV range. 
On low scales the effect of massive neutrinos is simply to 
lower the fluctuation level roughly as [E6| 



AP 



a. 



(12) 



0.020, Ho = 70 kms^Mpc" 



1.0, and t = 0, i.e. 



i.e. it is scale independent and therefore indistinguishable 
from changing the bias, b. This degeneracy can be bro- 
ken by precision measurements around the free-streaming 
scale where the break in the power spectrum occurs. A 
good example of how the mass limit on neutrinos can be 
tightened if bias is fixed comes from Ref. |50) . Here the 
mass limit comes from comparing the overall normaliza- 
tion of the spectra at COBE scales |3£| with those on 
cluster scales j57|. However, we believe that, at present, 
keeping bias as a free parameter yields a much more ro- 
bust constraint. To a much lesser extent the neutrino 
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TABLE I: The different priors on parameters other than Q v h 2 used in the analysis. 



prior type 


S2 m 


n b h 2 


h 


n 


r 


Q 


b 


CMB + LSS 


O.f-f 


0.008 - 0.040 


0.4- f.O 


0.66-f.34 


0-f 


free 


free 


CMB + LSS + BBN + Ho 


o.f-i 


0.020 ± 0.002 


0.70 ±0.07 


0.66-f.34 


0-1 


free 


free 


CMB + LSS + BBN + H + SNla 


0.28±0.f4 


0.020 ± 0.002 


0.70 ± 0.07 


0.66-f.34 


0-1 


free 


free 



mass is also degenerate with the Hubble parameter. 

It should be noted that there is little degeneracy with 
n s , the spectral index. In Refs. |2(| Q a significant 
degeneracy between fl v and n s was found when only LSS 
data is considered. However, this is broken when CMB 
data is included (as is also noted in Ref. |2^]), the reason 
being that changing n s affects both CMB and matter 
power spectra, not just the matter spectrum. 

Clearly, it would be desirable to fix the parameters 
with which the neutrino mass is most degenerate, fl m , 
b, and H n . As for fl m one can use the SNla result 
fl m = 0.28 ±0.14 which applies to a flat universe [Eg) . 
However, this value is not much more restrictive than 
what is found from the CMB±LSS data alone. Fixing 
the bias is much more difficult since it is not a physically 
well understood parameter. In Elgaroy et al. bias 
was kept as a free parameter, and we follow this line. Hq 
has been determined precisely by the HST Hubble key 
project to be H = 70 ± 7 kmV 1 Mpc" 1 |§. 

If all the above parameters are included in the Fisher 
matrix analysis the estimated la precision on m u is 1.8 
eV, equivalent to a 95% confidence limit of 3.6 eV. 



tations, roughly within the 68% confidence interval. 

We identify the 95% confidence limit on m v with the 
point where Ax 2 = 4. These limits are shown in Table 
II. For the most restrictive prior we find a 95% confidence 
upper limit of J2 m f — 2-47 eV. This is compatible with 
the findings of Ref. |l) who find that £ m v < 1.8 - 2.2 
eV for a slightly more restrictive parameter space. 

Based on the present analysis we consider Yl m v < 
3 eV (95% conf.) a robust upper limit on the sum of 
the neutrino masses. This corresponds roughly to the 
value found for the CMB±LSS data alone without any 
additional priors. Even though this value is significantly 
higher than what is quoted in Ref. |3lj|, it is still much 
more restrictive than the value m u < 4.4 eV |3£j found 
from CMB and PSCz |6^] data. As is also discussed 
in Ref. |3l| the main reason for the improvement is the 
much greater precision of the 2dF survey, compared to 
the PSCz data §§. 



III. NUMERICAL RESULTS 



The Fisher matrix analysis can only give a general idea 
of the constraints which can be found from a given data 
set. In reality the likelihood is non-Gaussian and away 
from the best fit point the formalism breaks down. In 
order to get reliable estimates it is necessary to perform 
a full numerical likelihood analysis over the space of cos- 
mological parameters. 

In this full numerical likelihood analysis we use a 
slightly restricted parameter space with the following free 
parameters: fl m , fib, Hq, n s , Q, b, and r. We restrict 
the analysis to flat models, flk = 0. This has very little 
effect on the analysis because there is little degeneracy 
between m v and fit: . In order to study the effect of the 
different priors we calculate three different cases, the pri- 
ors for which can be seen in Table I. The BBN prior on 
fib h 2 comes from Ref. [[30|. The actual maginalization 
over parameters other than fl u h 2 was performed using a 
simulated annealing procedure 

Fig. 2 shows x 2 for the three different cases as a func- 
tion of the m v . The best fit \ 2 values are 24.81, 25.66, 
and 25.71 for the three different priors respectively. In 
comparison the number of degrees of freedom are 34, 35, 
and 36, meaning that the fits are compatible with expec- 
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TABLE II: Best fit \ 2 and upper limits on Y ra„, ma x for the three different priors. 



prior type 


best fit x 2 


X) m ^max (eV) (95% conf.) 


CMB + LSS 


24.81 


2.96 


CMB + LSS + BBN + H 


25.66 


2.65 


CMB + LSS + BBN + H + SNIa 


25.71 


2.47 



IV. DISCUSSION 

We have studied cosmological constraints on the neu- 
trino masses from present CMB and LSS data. Initially 
a Fisher matrix analysis was performed which illustrates 
the main degeneracies of m v with other cosmological pa- 
rameters used in the analysis. From this simplified anal- 
ysis it was estimated that the precision on ^ m v should 
be roughly 3.6 eV at 95% confidence from CMB+LSS 
data alone. 

However, in order to obtain reliable estimates we per- 
formed a full numerical likelihood analysis. Using reason- 
able priors on fl m , Qbh and Hq we obtained a limit of 
Y^, m v < 2.47 eV at 95% confidence, while using no pri- 
ors on the CMB+LSS data yielded m v < 3 eV, again 
at 95% confidence. We believe this to be a robust upper 
limit. 

Our analysis shows, not surprisingly, that priors are ex- 
tremely important for parameter estimation of m v . Our 
most restrictive prior yields a result similar to that found 
by Ref . J3l| , while our no-prior case yields a significantly 
looser constraint. 

The Fisher matrix analysis showed that the parame- 



ter most degenerate with m v is the bias parameter, b. 
In order to obtain much stronger limits one must either 
determine b precisely in an independent way or obtain 
better LSS power spectrum statistics around the scale 
corresponding to the free-streaming scale for neutrinos, 
k ~ 0.02 - 0.03/i Mpc" 1 . The Sloan Digital Sky Survey 
will measure the power spectrum shape with higher 
precision on the relevant scale, and this data, combined 
with CMB data from the MAP experiment Q, will ei- 
ther push the limit by a factor of at least a few or indeed 
detect a non-zero neutrino mass directly. It was esti- 
mated in Ref. ]2(J that ^ m„ < 0.65 eV can be reached. 

Finally, we note that the present cosmological limit is 
significantly stronger than current laboratory limits. The 
Mainz tritium experiment |^2f currently quotes a 95% up- 
per limit to the v e mass of 2.2 eV, which translates to a 
sum of roughly 6.5-7 eV for the three mass eigenstates. 
As is also noted in Elgaroy et al. |3l), the cosmologi- 
cal limit is compatible with the very controversial detec- 
tion of neutrinoless double beta decay by the Heidelberg- 
Moscow experiment ]S5[ . If this finding is confirmed it 
would imply a sum of masses of order 1 eV, within range 
of the MAP+SDSS data. 
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